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Abstract 

The observed dark energy phenomenon is attributed to the presence 
of the zero-point fluctuations of matter fields. We show that due to the 
presence of the gas of virtual wormholes the zero-point energy is finite and 
forms the finite (of the Planckian order) value. The observed value of the 
cosmological constant is somewhat reduced due to the two effects. First is 
the renormalization of the gravitational constant which presumably forms 
the initial (local) value H? n f predicted by the inflationary scenario. And 
second is an additional reduction due to the presence of the gas of actual 
wormholes required by the dark matter phenomenon. We show that the 
Starobinsky model of inflation explains both, the inflationary stage in the 
past and the present day acceleration of the Universe. It also represents 
the so-called model of the eternal Universe. 

1 Introduction 

As is well known modern astrophysics (and, even more generally, theoretical 
physics) faces two key problems. Those are the nature of dark matter and dark 
energy. Recall that more than 90% of matter of the Universe has a non-baryonic 
dark (to say, mysterious) form, while lab experiments still show no evidence for 
the existence of such matter. Both dark components are intrinsically incorpo- 
rated in the most successful ACDM (Lambda cold dark matter) model which 
reproduces correctly properties of the Universe at very large scales (e.g., see [T] 
and references therein). The only failure of ACDM, i.e., the presence of cusps 
{rdm ~ l/ r ) i n centers of galaxies [2], can be cured if, instead of non-baryon 
particles, we shall use wormholes (see for detail [3]). Indeed, wormholes repre- 
sent extremely heavy (in comparision to particles) objects which at very large 
scales behave exactly like non-baryon cold particles, while at smaller scales (in 
galaxies) they strongly interact with baryons and form the observed [1] cored 
{pDM ~ const) distribution. It worth expecting that wormholes will play the 
central role in the explanation of the dark matter phenomenon. 

Save the dark matter component ACDM requires the presence ( ~ 70%) of 
dark energy (of the cosmological constant). Moreover, there are strong evidences 
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for the start of the acceleration phase in the evolution of the Universe [5]. In 
the present paper we shall try to demonstrate that wormholes play here the 
central role as well. It turns out that virtual wormholes are responsible for the 
formation of the renormalized (finite) value of the cosmological constant, while 
the observed value is somewhat reduced by the presence of the gas of actual 
wormholes. 

It seems that the observed small acceleration of the Universe is the start of 
a sequent inflationary stage predicted by the Starobinsky model [5] • Recall that 
the first inflationary model suggested in Ref. [B] is based on one-loop vacuum 
corrections to the Einstein equations 

R iw ~ \g^R = 7 {Tpv) (1) 

where the expectation value (T M1/ ) describes the energy of zero-point fluctuations 
of matter fields. In the absence of particles (vacuum state) 7 (T^) ~ A = const 
which immediately launches the inflationary phase. Such a phase was shown 
to be unstable in both directions. During the expansion the inflation ends due 
to the particle production, while in the backward direction the presence of an 
arbitrary small amount of matter (particles) destroys the inflationary phase as 
well (for on the last stage of the collapse the ordinary matter always prevails 
the cosmological constant). 

Wormholes add to the above scenario only two features. First, virtual worm- 
holes remove divergencies and make the energy of zero-pint fluctuations be finite. 
This, in principle, allows to account for higher order corrections to the above 
equation. And the second feature is that wormholes always remove some por- 
tion of degrees of freedom and, therefore, they can be described by ghost fields 
(by the analogy with Fadeev-Popov ghosts). The inflationary phase is always 
accompanied with the production of wormholes (or ghost particles) as well, 
which leads to an additional reduction (i.e., a partial decay) of the cosmological 
constant (zero-point energy). On the subsequent stage wormholes created can 
merge reheating thus the rest matter. This gives the standard Friedman stage 
of the evolution which lasts until the cosmological constant prevails again the 
matter (as we do observe now) and all the picture repeats. It is impossible to 
say on which cycle of this eternal process we are. Thus the Starobinsky model 
seems to give the simplest model of the eternal Universe. 

Since the divergencies in QFT (quantum field theory) come from the ultra- 
violet behavior which involves extremely small scales, we may expect that effects 
of general relativity are not essential here and in the present paper we merely 
neglect them to the leading order. We also use everywhere the Planckian units. 
In the present paper we also use only mass- less wormholes [HI 13]. In general 
massive wormholes may also appear but we expect that they are essentially 
suppressed (from the energetic standpoint). Moreover, we may expect that 
basic topological effects will not change. We point out that wormhole rest 
masses do appear in the presence of particles created as it was explained in [3J 
which leads to some additional non-linear phenomena on the inflationary stage. 
This however requires an independent investigation. 
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The paper is organized as follows. In Sec. 2 we present the construction of 
the generating functional with virtual wormholes taken into account. In Sec. 3 
we investigate properties of the two-point Green function. We show that the 
presence of the gas of virtual wormholes can be described by the topological bias 
[7] exactly as it happens in the presence of actual wormholes [HI |3] . We demon- 
strate that the mean value for the bias defines the cutoff function in the space 
of modes. In Sec. 4 we explicitly demonstrate that for a particular set of worm- 
holes the bias defines not more than the projection operator on the subspace of 
functions obeying to the proper boundary conditions at wormholc throats. The 
projective nature of the bias means that wormholes merely cut some portion of 
degrees of freedom (modes) . Phenomenologically it means that wormholes can 
be described by the presence of ghost fields which compensate the extra (cut 
by wormholes) modes. In Sec. 5 we show how the cutoff expresses via some 
dynamic parameters of wormholes. The exact definition of such parameters we 
leave for the future investigation. In Sec. 6 we demonstrate that wormholes lead 
to a finite (of the planckian order) value of (T^) which requires considering the 
contribution from the smaller and smaller wormholes with divergent density 
n —> oo. In Sec. 7 we define the renormalization of the gravitational constant 
which leads to a reduction of the observed value of A. Recall that inflationary 
scenarios require \/A/3 = Hi n f ~ ~ AT/T which in our approach means 
that 7 ~ 10 _1 °7o, where 70 is the naked gravitational constant. In Sec. 8 we 
consider actual wormholes and show that they produce a negative contribution 
to the cosmological constant. Thus, the creation of wormholes during the in- 
flationary stage leads to a partial decay of the initial cosmological constant. In 
Sec. 9 we repeat basic results an discuss some perspectives. 

2 Generating function 

The basic aim of this section is to construct the generating functional which can 
be used to get all possible correlation functions. Consider the partition function 
which includes the sum over topologies and the sum over field configurations 



For the sake of simplicity we shall use from the very beginning the Euclidean 
approach (e.g., see Refs. [S]-[H] and references therein), i.e., the action in the 
form 



and the notions (Jtp) = J J (x) tp (x) d x, where the integral is taken over phys- 
ically admissible portion of i? 4 only. From the variation principle SS = we 
may find the equation of motions for the field tp 




(2) 




(3) 



Atp + J = 



(4) 
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which has the solution (A 1 is the Green function obeying to a proper boundary 
conditions) 

ip = A- 1 J. (5) 

2.1 Sum over fields 

We fix the topology of space by placing a set of wormholes with parameters 
£i (i.e., = (di, , R^) , where a is the throat radius and R^ are positions 
of throats in space). For general properties of a wormhole see Ref.[3]. Then 
we consider the sum over field configurations ip, which can be replaced by the 
integral 

Z*(J)= f [LV]eH^M^). (6) 



Upon the simple transformations 

\ (ipAip) (J<p) = \ (ipA?) - \ (ji- 1 j) , (7) 
where tp = ip — A J, we cast the partition function to the form 

Z* = [ [D0\ e^ A ^( jA ~ 1J ) = Z (A)e-H jA - lj ), (8) 



where Zq(A) = J [Dip] e 2 ( v ^ v ) is the standard expression and A 1 = A 1 (£i,...,£iv) 
is the Green function for a fixed topology, i.e., for a fixed set of wormholes 
£ij — >&v ■ 

In the case of a fixed topology the generating functional ([5])-© allows us to 
construct the perturbation scheme, when we add to the action ([3]) the interaction 
term AS (ip) , by means of using the obvious expression 

Z (J) = e^i^Z* (J) = £ ( A s (-^))" Z* (J) (9) 

and generate all possible momenta (the higher order Green functions G s (xi, X2, ...£ s )) 
as 



which depend on parameters of wormholes, i.e., G s = G s (£i, £n). 



(10) 



2.2 Sum over wormholes 

Consider now the sum over topologies r. To this end we restrict with the sum 
over the number of wormholes and integrals over parameters of wormholes: 



where in general the integration over parameters is not free (e.g., it obeys the 



obvious restriction 



> 2<Xj). This defines the generating function as 



Z to tai (■/) = £/ (II^ ) Z Q {A)e-^ jA ' Xj ). (12) 

N J \i=l J 

The sum over topologies assumes an additional averaging out for the Green 
functions (fT0|) with the measure djiM = p (£, N) d N £, where 



(0) 



J total 



(13) 



which obey the obvious normalization condition YIn I ^MJV = YIn P n = 
The averaging out over topologies assumes the two stages. First we fix the total 
number of wormholes N and average over the parameters of wormholes £ (i.e., 
over parameters of a static gas of wormholes in i? 4 ). Then we sum over the 
number of wormholes N (the so-called big canonical ensemble) . 

The basic difficulty of the standard field theory is that the perturbation 
scheme based on the decomposition © leads to divergent expressions. This 
remains true for every particular topology of space (i.e., for any particular set 
of wormholes), since there always exists a scale below which the space looks like 
the ordinary Euclidean space. What we expect is that the sum over all possible 
topologies will remove such a difficulty. 



3 The two-point Green function 

From (|5J), ©, and (JTDJ) we see that the very basic role in QFT plays the two 
point Green function. Such a Green function can be found from the equation 

AG(x,x')=-5(x-x') (14) 

with proper boundary conditions at wormholes, which gives G = A^ 1 . Now let 
us introduce the bias function N(x, x') as 

G(x,y) = f G (x,x')N(x',y)dx', (15) 

where Gq(x,x') is the ballistic (or the standard Euclidean Green function) and 
the bias can be presented as 

N(x,x') = 5{x-x') + ^b l 5{x-x i ) (16) 

i 

where 6j are fictitious sources at positions x\ which should be added to obey 
the proper boundary conditions. 
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In the simplest approximation (of a rarefied gas) the bias function can be 
expressed via the bias for a single wormhole. Indeed if we use the decomposition 
N(x, x') — 5(x, x') + b (x, x'), then 



b(x,x')=Nj h(x,x',OF(Odt (17) 

where bi(x,x',^) is the bias function b for a single wormhole, N is the total 
number of wormholcs, while 

1 N 

i=l 

is the density of wormholes in the configuration space £. In what follows we 
shall use the notion 

ni (x, x',£) = S(x - x 1 ) + Nbi (x, x',£) 
which, in virtue of J F(^)d£_ — 1, defines the total bias function as 

N(x,x') = f ni (x,J,S)F(Qdt. (19) 



We point out that the superposition of wormholcs in (|19[) represents only an 
approximation (which works when we retain only the first order images by the 
parameter <C 1 (e.g., see Ref. [5]), while in the most general case the 

superposition does not work and N(x, x') = N(x, a/, £i, ...£/v))- 

Now taking into account that due to <jT5j) the total Green function can be 
expressed via the auxiliary Green function G\ (x, y, £) — J Gq(x, x')ni(x\ y)dx' 
(which in the case N = 1 gives the Green function for a single wormhole) as 
G = J G\ (£) F(£)d£. Let us introduce the notion 

[jGj] = J (jd!j) F($dZ, (20) 

then for the sum over wormholes (|11[) we may replace 

Ztotai ( J) = J [DF] Z {G- l )e-?{ j6j } (21) 

N 

where / [DF] =J2I U d ^- 

N i=l 

The action (J3]) remains invariant under translations x! = x + c with an 
arbitrary c which means that the measure (|13[) does not actually depend on the 
position of the center of mass of the gas of wormholes and, therefore, we may 
restrict ourself with homogeneous distributions F(£) of wormholes in space only. 
Indeed, we may define d N £ — d N £'d 4 c, while the integration over d 4 c gives the 
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volume of R 4 i.e., / d 4 c = L 4 = V which disappears from (fT3f due to the 
denominatofl In what follows we shall omit the prime from 

Let us consider the Fourier representation N(x, x', £) —> N(k, k' , £) which in 
the case of a homogeneous distribution of wormholes gives N(k, k') = N(k, £,)S(k— 
k'), then we find 

G(k) = G (k)N(k,$ 



and the Green function can be taken as 



k-Z J- rr,i 



Then for the total partition function we find 

Z to tal(J)= I [DN{k)]e~^ N ^Z (N(k))e~ h ^^^ V,kf ) d \ (23) 



where [DN] — J\dNk and <r(N) comes from the integration measure (i.e., from 

k 

the Jacobian of transformation from F (£) to N (k)) 

g tr(iV) f [ DF ]S(N(k)-N(k,t)). 



We point out that a(N) can be changed by means of adding to the action ([3]) 
of an arbitrary " non-dynamical" constant term which depends only on topology 
(wormholes) S -> S + AS(N (jfe)). The multiplier Z (N) defines the simplest 
measure for topologies which is given by 

Z = cxp (- f Ad 4 x] = exp J ±i-^-r / In ( , 5 - 9 ) N{k)d 4 k \ , (24) 



"2(27r) 4 y \k 

where ± stands for Bose/Fermi statistics of the field ip. Now by means of using 
the expression (|23[) and (|10p we find the two-point Green function in the form 

G ( fc ) = T^T ( 25 ) 
k + ra l 

where N(k) is the cutoff function (the mean bias) which is given by 

N(k) =J2 [ dmN (k, = 7 1 M I [DN] e-^Z (N) N (k) . 

n J * total IUJ J 



N 



At the present stage we still cannot evaluate the exact form for the cutoff 
function N(k) in virtue of the ambiguity of AS(N (k)) pointed out. Such a term 
may include two parts. First part AiS describes the proper dynamics of worm- 
holes and should be considered separately, while the second part may describe 



1 Technically, we may first restrict a portion of Ft 4 in the to a finite volume V and 
then in final expressions consider the limit V — > oo (which represents the standard tool in 
thermodynamics and QFT). 
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"external conditions" for the intensity of topology fluctuations. Actually the 
last term can be used to prescribe an arbitrary particular value for the cutoff 
function N(k) = f(k). Indeed, the "external conditions" can be accounted for 
by adding the term A 2 S = (X,N) = f A (fc) N (k) d 4 k, where A (k) plays the 
role of a specific chemical potential which implicitly depends on / (k) through 
the equation 

f(k) = -= * / [DN]e^-^N(k), 

A total (A UJ J 



where the value of Zo (|24|) is included into the value of the chemical potential 
A (fc). From the QFT standpoint such a term leads merely to a renormalization 
of the cosmological constant. By other words the intensity of topology fluctu- 
ations (i.e., the cutoff function) is driven by the cosmological constant A and 
vise versa. 



4 Topological bias as a projection operator 

By the construction the topological bias N (x, x') plays the role of a projection 
operator onto the space of functions (a subspace of functions on i? 4 ) which 
obey the proper boundary conditions at throats of wormholes. This means 
that for any particular topology (for a set of wormholes) there exists the basis 
{fi ( x )} m which it takes the diagonal form N (x, x') = ^ifi ( x ) ft ( x ') with 
eigenvalues Ni — 0,1 (since N? = Ni). In this section we illustrate this simple 
fact (which is probably not obvious for readers) by the explicit construction of 
the reference system for a single wormhole when physical functions become (due 
to the boundary conditions) periodic functions of one of coordinates. 

Indeed, consider a single wormhole with parameters £ (i.e., £ = (a, R + , R~), 
where a is the throat radius and are positions of throats in spacqj) . Consider 
now a particular solution (f> to the equation A0 O = (harmonic function) for 
i? 4 in the presence of the wormhol^l, which corresponds to the situation when 
throats possess a unit charge/mass but those have the opposite signs. Now define 
the family of lines of force x (s, xq) which obey the equation ^§ = — V0o ( x ) with 
initial conditions x(0) = xq. Physically, such lines correspond to lines of force for 
a two charged particles in positions R^ with charges ±1. We note that all points 
which lay on the trajectory x (s, xq) may be taken as initial conditions and they 
define the same line of force with the obvious redefinition s — > s — sq . By other 
words we may take as a new coordinates the parameter s and portion of the 
coordinates orthogonal to the family of lines Xq . Coordinates Xq can be taken 

2 In general, there exists an additional parameter Ug which defines a rotation of one of 
throats before gluing. However, it does not change the subsequent construction. There always 
exists a diffeomorfic map of coordinates x' = h(x) which sets such a matrix to unity. 

3 Instead of the construction used here one may use also another method. Indeed, consider 
two point charges, then the function </>o = 1/ (x — x+) 2 — l/(x — X-) 2 can be taken as a new 
coordinate. Wormhole appears when we identify (glue) surfaces ipo = iw. We point out that 
such surfaces are not spheres, though they reduce to spheres in the limit \x+ — X- —> oo or 
U) — > oo. 
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as laying in the hyperplane R 3 which is orthogonal to the vector d = R~ — R + 
and goes through the point X = (R~ + R + )/2. Let s^ 1 (xq) be the values of the 
parameter s at which the line intersects the throats i? . Then instead of s we 
may consider a new parameter 9 as s (9) = s~ + (s + — s~) 9/2tt, so that when 
9 = 0,2n the parameter s takes the values s = s~ , s + respectively. The gluing 
procedure at throats means merely that we identify points at = and 9 = 2n 
and all physical functions in the space i? 4 with a single wormhole £ become 
periodic functions of 9. Thus, the coordinate transformation x = x(9,Xq) gives 
the map of the above space onto the cylinder with a specific metric dl 2 = 
(dx(9,XQ^ = g a /3dy a dyP (where y — (6, Xq)) whose components are also 
periodic in terms of 9. Now we can continue the coordinates to the whole space 
R 4 (to construct a cover of the fundamental region e [0, 27r]) simply admitting 
all values — oo < 9 < +oo this, however, requires to introduce the bias 

+ oo 

— 5(9-9')^ N {9 -9')= (9 - 6' + 2nn) , 

V " n=-oo V " 

since every point and every source in the fundamental region acquires a count- 
able set of images in the non-physical region (inside of wormhole throats). Con- 
sidering now the Fourier transforms for 9 we find 



-f-oo 

N(k,k')= Hk-n)S(k-k'). 



We point out that the above bias gives the unit operator in the space of periodic 
functions of 9. From the standpoint of all possible functions on R A it represents 
the projection operator ./V 2 = N (£) (taking an arbitrary function / we find 
that upon the projection /jy = Nf Jn becomes a periodic function of 9, i.e., only 
periodic functions survive). 

The above construction can be easily generalized to the presence of a set 
of wormholes. In the approximation of a dilute gas of wormholes we may ne- 
glect the influence of wormholes on each other (at least there always exists a 
sufficiently smooth map which transforms the family of lines of force for "in- 
dependent" wormholes onto the actual lines). Then the total bias (projection) 
may be considered as the product 



N tota i{x,x')= / \ X\y/gld?yi ) N {^,x, yi ) N (^, yi ,y 2 ) ...N ^ N ,y N _ u x') 




where N(£i,x,x') is the bias for a single wormhole with parameters Every 
such a particular bias N (£j, x, x') realizes projection on a subspace of functions 
which are periodic with respect to a particular coordinate 9i (x) , while the to- 
tal bias gives the projection onto the intersection of such particular subspaces 
(functions which are periodic with respect to every parameter #,). 
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5 Cutoff 



The projective nature of the bias operator N (x, x') allows us to express the 
cutoff function N(k) via dynamical parameters of wormholes. Indeed, consider 
a box L in i? 4 and periodic boundary conditions which gives k = 2im/L (in 
final expressions we consider the limit L — > oo, which gives ^2 k — > j^js J d 4 k). 
And let us consider the decomposition for the integration measure in (|23p as 

a (N) = <t q + Ai (k) N (k) + i A 2 ( fc > k ') N ( fc ) N i k ') + - 

where Ai (k) includes also the contribution from Zq (k). We point out that this 
measure plays the role of the action for the bias N (k). Indeed, the variation of 
the above expression gives the equation of motions for the bias in the form 

J2^2 (k,k')N(k r ) = -X 1 (k) 

k' 

which can be found by considering the proper dynamics of wormholes. We 
however do not consider the problem of the dynamical description of wormholes 
here and leave this for the future research. Then in the first approximation we 
may retain the linear term only. Then taken into account that N (k) =0,1 
(N 2 = N) we find 

(26) 

The simplest choice gives merely Ai (fc) = — In Zq (fc), where the sum is taken 
over the number of fields and Zq (k) is given by Zq (k) = ^ix/ (k 2 + to 2 ). In 
the case of a set of massless fields we find N(k) = Z (k) / (1 + Z (k)) where 
Z (k) = (^/tt /k) and a is the number of degrees of freedom. To ensure the 
absence of divergencies one has to consider the number of fields a > 4 [7]. 
However, as we show in the next section such a choice cannot be correct. Indeed, 
while its behavior at very small scales (i.e., when exceeding the Plankian scales 
Z (k) <1 and N(k) = Z (fc)) may be physically accepted, since it produces some 
kind of a cutoff, on the mass-shell k 2 + m 2 — > it gives the behavior N(k) — > 1 
which is merely incorrect (the true behavior is N(k) —> const < 1). 



6 Cosmological constant 

In this section we consider the renormalization of the stress energy tensor. Let 
N be the fixed number of wormholes. Then the stress energy tensor can be 
obtained directly from the two-point green function (|2"2"1) . (|2~5j) as 

- (T afj (a:)) = Um (d a d' p - l -g afi (d% + to 2 )") (G (x, x', £)) . (27) 
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By means of using the Fourier transform G [x, x', — / e ik ^ x x ^ s , £) (f^yr 
and the expression f|22[) . ([25]) we arrive at 

-(Ta fi {x))=-\g afi J (l + ^-^{NM)^, (28) 

where the property J k a kpf (fc) d 4 fc = |g a/ 3 / k 2 f (k) d A k has been used. 

For the sake of simplicity we consider the massless case. Then by the use of 
the cutoff N{k) — 7r Q / 2 / (tt"/ 2 + fc Q ) from the previous section we get the finite 
value (a > 4 is the number of the field helicity states) 

1 f W 2 d 4 k 1 /a-4\ /4\ , , 

A = W -Y2 V {—) r U) • ^ 

Since the leading contribution comes here from very small scales, we may hope 
that this value will not essentially change if the true cutoff function changes the 
behavior on the mass-shell as k — > (e.g., if we take Ai (k) = — ^lnZ (k) + 
S\ (k) with 5X (k) < In Z (k) as k > 1). 

To understand how wormholes remove divergencies, it will be convenient 
to split the bias function into two parts N (k,£) = 1 + b(k, £), where 1 corre- 
sponds to the standard Euclidean contribution, while b (k, £) is the contribution 
of wormholes. The first part gives the well-known divergent contribution of 
vacuum field fluctuations (t^\ — Ag a p with A — > +oo , while the second part 

remains finite for any finite number of wormholes and, due to the projective 
nature of the bias described in the previous section, it partially compensates 
(reduces) the value of the cosmological constant, i.e., (AT Q( g) = —5Ag a p, where 
5A = J2n PnSA (N) and SA (N) is a finite contribution of a finite set of worm- 
holes. 

6.1 The long- wave behavior of the bias (rarefied gas ap- 
proximation) 

To illustrate the compensational role of virtual wormholes we consider now the 
bias for a particular set of wormholes. For the sake of simplicity we consider 
the case when m = 0. Consider the Green function for the Laplace equation 

-AG {x,x')=6(x-x') 

in the presence of a single wormhole. A single wormhole can be viewed as a 

couple of conjugated spheres S± of the radius a with a distance d = R+ — R- 

between centers of spheres. So the parameters of the wormhole are0 £ = 
(a, -R-i-, i?_). The interior of the spheres is removed and surfaces are glued 
together. The Green function for the Euclidean space is merely Gq(x,x') — 



4 We do not account here for the additional parameter (rotation matrix U), since in the 
approximation used such parameter gives not essential contribution. 
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— ? 7^j~7p (Co (fc) = l/^ 2 )- The actual Green function can be found by means 
of using the image method (the straightforward generalization of results in Rcf. 
[8]). Then the proper boundary conditions (the actual topology) can be ac- 
counted for by adding the bias of the source 



6(x - x') ->• S(x - x') +b(x,x') 



(30) 



In the approximation a/d <C 1 (e.g., see for details Ref.[8]) the bias takes the 
form 



h = a 2 



1 



1 



(R- - x'f (R+ - x'f 



5(x- R + )- S(x- R_) . (31) 



This form for the bias is convenient when constructing the true Green function 
and considering the long- wave limit, however it is not acceptable in considering 
the short-wave behavior and vacuum polarization effects. Indeed, the positions 
of additional sources are in the physically non-admissible region of space (the 
interior of spheres S±). To account for the finite value of the throat size we 
should replace in (|3"Tj) the point-like source with the surface density (induced on 
the throat) i.e., 

1 



6(x-R±) 



x — R± 



a). 



(32) 



2ir 2 a 3 

Such a replacement does not change the value of the true Green function, how- 
ever now all extra sources are in the physically admissible region of space. 

We see that the bias function for the gas of wormholes in this approximation 
is additive, i.e., 



hiatal (x,x')=Y, bl {X,X',^)=N J b X (x, x' , 0^(£R- 



(33) 



For a homogeneous and isotropic distribution F(£) — F (a, d), then for the bias 
we find 



btotai (x - x 1 ) 



1 



7r 2 a V R 



1 



1 



X — x' — i?4 



-a)NF {a,\R+-R-\)d£ 
(34) 



Consider the Fourier transform F (a, X) = j F {a.k) e tkx -p^yr and using the 
integral ^ = J ^e~ lkx d4fc 



(2tt) 



btotai (k) = 2N 



^4we find for b (k) — J b (x) e lkx d A x the expression 



,4tt 2 



(F(a,k)-F(a,0)) 



Ji (ka) 
ka/2 



da. 



(35) 



Consider now a particular form for F(a,X), e.g., 



NF (a,X) 



2tt 2 ? 



r <5 (a - a ) S (X - r ) 



(36) 



where n = N/V is the density of wormholes. In the case N — 1 this func- 
tion corresponds to a single wormhole with the throat size ao and the distance 
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between throats ro = \R+ — R-\. We recall that the action ^ remains in- 
variant under translations and rotations which straightforwardly leads to the 
above function. Then NF (a, k) = J NF (a, X) e %kx d A x reduces to NF (a, k) = 

(a - a ). Thus from K3SJ we find 

w k 2 V kr /2 J ka/2 v ' 

As fc -> we get J x (fcr ) /^f 1 « 1-| (^l* 1 ) +••• which gives 6 (fc) « -7r 2 na 2 r 2 + 
.... In a more general case we find that on the mass-shell {k — > 0) b(k) s» 
— J 7r 2 a 2 rgn (a, ro) dadr$ + where n (a, ro) is the density of wormholes with 
a particular values of a and ro, and for the cutoff function (|26j) we get 

_ r r 2 

AT(fc) -> 1 - / vrVn (a, r ) ^dadr < 1. (38) 

Thus, we see that on the mass-shell the presence of virtual wormholes di- 
minishes the value of the cutoff function, which should lead to a specific renor- 
malization of the field, rest mass, and charge values. 

6.2 Renormalization of the cosmological constant 

Consider now the contribution of virtual wormholes to the cosmological constant 
which is given by 5A(N) = -|J"&W(Qt = ~k b total (0). Then from the 
expressions (|34| and (|36|) we get 

&t oio ; (0) = - 27r4 ^3 / (i - |r) - «)*(l#f --R-I - r )d 4 R-d 4 R+ 



which gives 



where 



</ \ J / a 2 sin 2 0d0 



r o 7 I" 7o fl2 + 2aro cos # + 



'o 

For a/ro <C 1 (we recall that by the construction a/ro < 1/2) this function has 
the value / ~ a2 / r o- Thus, for the contribution of wormholes we find 

If, ^ 1 _ f f^_^ ■ ■ n 



«5A = -y n(o,r ) I 1 - / I -J j dadr = - (1 - (/)) . 

From the above expression we see that to get the finite value of the cos- 
mological constant Ao = A — 6 A < oo one should consider the limit n — > oo 
(infinite density of virtual wormholes) which requires considering the smaller 
and smaller wormholes. From the other hand we have the obvious restriction 
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f 2n^j" a 4 dadro < 1, where ^-a 4 is the volume of one throat (wormholes cannot 
cut more, than the volume of space). Therefore the infinite density of worm- 
holes remains consistent with a finite renormalization (|38p (e.g., for the density 
of wormholes n(a) it is sufficient to have the behavior n(a) ~ 1/a as a — » 0). 
We also point out that the problem to determine the exact value for the cosmo- 
logical constant Ao requires to involve the simultaneous renormalization of all 
other fundamental constants (rest masses, charge values, and the gravitational 
constant). 



7 Renormalization of charge values 

The absolute value of the renormalized cosmological constant (e.g., see (p25j) ) 
does not give the observed value yet. Indeed the behavior of the mean cutoff 
function at the mass-shall N(k) ([55)1 defines some renormalization (screening) of 
the gravitational constant (of charge values) as well, while (T a p (x)) represents 
the source in (fT]). This leads to a somewhat reduced value of the cosmological 
constant. To find the observed value we have to account for two points. First is 
that the gas of virtual wormholes is actually dense (n — > oo). And the second is 
that the polarization at wormholes has some part which leads to anti-screening 
(for the sake of simplicity we have neglected this part in (l3"Tj) e.g., see for detail 
Ref.g]). 



7.1 dense gas 

The expression (|38p was obtained in the approximation of a rarefied gas. Indeed, 
only in this case the total bias obeys the superposition property (|33|) . The dense 
gas can be accounted for as follows. Consider the equivalent description which 
is the introduction of the topological permeability e, i.e. the modification of 
the equation in the form AsG(x, x') = — An 2 8(x — x'). In the case of the 
homogeneous and isotropic gas the relation between the bias b and e is simple in 
the Fourier representation, i.e. e(k) = 1/(1 + 6 (fc)) s» 1 — b(k). In the situation 
when e = const (e.g., at very large scales k — > 0) the topological permeability 
renormalizes merely the value of a source AG(x, x') — (4ir 2 /e)6(x — x') (or 
equivalently the value of the interaction constant 7 — > 7/e ). Then (for a 
constant e) the permeability of a dense gas can be obtained in the standard 
way. Indeed, let present e = 1 + 47r 2 %, where \ is the topological susceptibility 
of space. In linear approximation we get 

f r 2 
47T 2 Xo ~ / 7r 2 a 4 n (a, r ) -^dadr , 



see (|38[) . Then for a dense gas it is related to the linear susceptibility xo as 
X = Xo/ (l ~ 7I ' 2 Xo), e.g., see part 4.6 in Ref. [T3]. Thus finally, we find the 
renormalization of the physical value of the interaction constant in the form 

7 (l-7T 2 Xo) 
Iphys = -= 7 -1 : o 2 < 7- 

£ 1 + 3tt 2 xo 
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7.2 renormalization of volume 



The anti-screening part of the bias somewhat diminishes the reduction of the 
cosmological constant. It can be accounted for by the renormalization of the 
physical volume [5] . Indeed, wormholes merely cut some portion of the coordi- 
nate volume which is given by 



Vphy S — V I 1 - J 7r a n (a, r ) dadr 

where the second term describes the volume of interiors of wormholes. Suppose 
that we work in the standard Euclidean space R 4 (we do not aware that the 
actual space has a smaller volume) . Then we should introduce the apparent (or 
observed) interaction constant as 

V 

^/apparent — ^fphys T7 ^ T 

*phys 

which accounts for the fact that the number density of strength lines remains 
the same, while the total number of strength lines fictitiously increases when we 
continue fields to the non-physical regional Thus for the observed value of the 
cosmological constant we get 

V 

Kbs = IphysT^ A «A~1. 

vphys 

8 Dark energy from actual wormholes 

Consider now the contribution to the dark energy from the gas of actual worm- 
holes. Unlike the virtual wormholes, actual wormholes do exist at all times and, 
therefore, a single wormhole can be viewed as a couple of conjugated cylinders 
T± = S± X R . So that the number of parameters of an actual wormhole is less 
•q = (a, r+, r_), where a is the radius of S± and r± £ R 3 is a spatial part of R±. 

For rigorous evaluation of dark energy in this case we, first, have to find the 
bias bi {x,x',rj) analogous to (l3Tj) for the topology i? 4 /(T+ UT 3 ). There are 
many papers treating different wormholes in this respect (e.g., see Ref. |15) and 
references therein). However, in the present paper for an estimation we shall 
use a more simple trick. 



5 Such an effect can be easily analyzed on the simplest example of a cylinder R 2 X S 1 . In 

terms of R? the number density of degrees of freedom for the cylinder is always N (k) < 1. 

Indeed, consider the finite volume V € R 3 and plane waves fu = —j= e~ lkx which means that 

vV 

k = 2nn/L (V = L 3 ). Then £k = / J^pjj^k = J J^ N ( k ) d3fc > i- e -> for r3 we have N = 1 - 

In the case of the cylinder we find = f 4 d 3 k, for small distances k 2> k* = l/£ (£ is 

2 

the radius of the cylinder) and y\ = / L ^ (fik as k* 2> k 2> ko, where fco = 2n/L. This 
defines N c (k) = V phys /V = i/L < 1 as k > k* and N c (k) = ir/ (Lk) < 1 as k* > k > fc - 
Thus if we continue the cylinder to the whole space R? we set 7V C = 1 as k 3> k* , while for 
scales k* 3> k we will get N c = it/ (£k) ^> 1 which may be interpreted as the presence of dark 
matter (the scale dependent renormalization of the apparent mass M a p P arent (R) ~ MoRk* 

ED- 
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8.1 Beads of virtual wormholes 

Indeed, instead of the cylinders T± we consider a couple of chains (beads of vir- 
tual wormholes T± — > Ll n S± „). Then the bias can be written straightforwardly 



h (x,x',ti) 



E 



i 



i 



i 



2tt 2 c 



(i?_,„ - x'Y (R 



x-R 



a)-S( 



R- 



+ ,n 



a) 



x'Y 



(39) 



where R±, n = {t n ,r±) with t n = to + 2£n and I > a is the step. We may 
expect that upon averaging over the position to € [— the bias for the beads 
will reproduce the bias for cylinders T± ( at least it looks like a very good 
approximation) . We point out that the averaging out ^ J_ g dto and the sum 

Y^n=-oo reduces to a single integral ^ dt of the zero term in (13T))) . And 
moreover, the resulting total bias corresponds merely to a specific choice of the 
distribution function F(£) in (|33[) . Namely, we may take 



t-)f(\r+-r-\,a), 



where R± = (t±, r±) and / (s, a) is the distribution of cylinders, which can be 
taken as (n is 3-dimensional density) 



n (a) 
47rrg 



ro). 



Using the normalization condition J NF(£)d£ — N we find the relation A = 
■^nV — nV, where n is a 4-dimensional density of wormholes and 1/ (2^) is 
the effective number of wormholes on the unit length of the cylinder (i.e., the 
frequency with which the virtual wormhole appears at the positions r± ). This 
frequency is uniquely fixed by the requirement that the volume which cuts the 
bead is equal to that which cuts the cylinder |7ra 3 = \o^^g (i.e., 21 = and 
n = ^—n). Thus, we can use directly expression (I3"5t and find (compare to (|37[) 



b{k) 



2n (a) a 



,4tH 



1 



sin |k| ro \ J\ (ka) 



|k|r 



ka/2 



-da, 



(40) 



where k = (fco,k). Here the first term merely coincides with that in (|37l) and, 
therefore, it gives the contribution to the cosmological constant A = — n/2 = 
—An/ (37ra), while the second term describes a correction which does not reduce 
to the cosmological constant and requires a separate consideration. 



8.2 Stress energy tensor 

From (|27|) we find that the stress energy tensor 



(AT Q/3 (x)) 



d A k 



(41) 
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reduces to the two functions 

Too = e = Ai — 2^ 

T 2i = p*y, p = ^A 2 - -A* 
where e + 3p = —fjt and Ax + A2 = /i and these functions are 

/-fco 2 ,^ > /-|k| 2 ,rf 4 fc [ h d*k 



(2tt)*' 7 (27rf 7 (2tt) 

By means of the use of the spherical coordinates k^/k 2 — cos 2 6, |k| 2 /fc 2 = 
sin 2 0, and d 4 fc = An sin 2 9k 3 dkd8 we get 

* - f 1 - 2 fe) 2 * (£)) ■ « - ^ f 1 - 1 fe) 2 * 



where 

sin 6dxd6. 



(a\ 2 f f°° . , . a , Jl (jj ^cos 2 
For - « 1 we find 

Thus, finally we find 

4n / 4/a\ 2 \ 4H /. 4 / a x jX 



2 37ra ' ^ 1 3 \ ro / y 37ra 1 3 \ ro 

Which upon the continuation to the Minkowski space gives the equation of state 
in the fornix 



We point out that on the contrary to the virtual wormholes the energy density 
of actual wormholes is always finite — e < 00. 



9 Summary 

Thus, we see that virtual wormholes indeed lead to the regularization of all 
divergencies in QFT and form the local value of the cosmological constant. 

6 An arbitrary gas of wormholes splits in fractions with a fixed a and ro. 
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This automatically distinguishes the Starobinsky model of inflation [5] as the 
model of the eternal Universe. It is impossible to say at which inflationary cycle 
we are. 

In the case of actual wormholes the leading contribution to the stress energy 
tensor gives the negative cosmological constant. This happens due to the general 
fact that wormholes always remove some portion of field modes. Phenomenolog- 
ically, such a removal can be prescribed to the presence of ghost fields. Virtual 
ghost particles reflect virtual wormholes, while real ghosts will describe worm- 
holes created. On the inflationary stage production of ghost particles (of actual 
wormholes) essentially reduces the physical value of the cosmological constant 
which leads to a somewhat more rapid escape from the inflationary regime. We 
also point out that particle production during the inflation automatically gen- 
erates huge rest masses of wormholes [3] which should lead to some complex 
picture for the escape from inflation. We also point out that the relating (see 
Sec. 6) dynamic parameters of wormholes (or ghost fields) to the cutoff function 
requires further investigation. 

We also point out that in the case of virtual wormholes the topological 
polarization effects (e.g., see the analogous consideration of actual wormholes in 
Ref. [3]) induce "dark charge" values on throats. This explains the well-known 
phenomenon of the dark charge (or the Higgs sector) in particle physics. We 
however consider this problem elsewhere. 
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